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ABSTRACT 

The problem of a cylindrical shell containing a circumferential 
through crack is considered by taking into account the effect of trans- 
verse shear deformations. The formulation is given for a specially 
orthotropic material within the confines of a linearized shallow shell 
theory. The particular theory used permits the consideration of all 
five boundary conditions regarding moment and stress resultants on the 
crack surface. Consequently, aside from multiplicative constants repre- 
senting the stress intensity factors, the membrane and bending components 
of the asymptotic stress fields near the crack tip are found to be iden- 
tical. Since these two fields are coupled in the shell problems, the 
difficulty in justifying the superposition of membrane and bending 
stresses and, particularly of the related stress intensity factors when 
applying the results to fracture problems has thus been removed. The 
stress intensity factors are calculated separately for a cylinder under 
a uniform membrane load, and that under a uniform bending moment. The 
stress intensity factors are calculated for values of 0 < X 2 = [12(l-v 2 )]' t 
a//Rh<10 and 0<a/h<10, where 2a, h, and R are, respectively, the crack 
length, the thickness, and the (mean) radius of the cylinder. Sample 
results showing the nature of the out-of-plane crack surface displace- 
ment and the effect of the Poisson's ratio are also presented. 

1. INTRODUCTION 

Depending on the mechanical behavior of the material and the nature 
of the loading and environmental conditions, in designing the line pipes 


and cylindrical containers it is often necessary to consider fatigue crack 
propagation and fracture among the possible modes of failure. This requires, 
in addition to the application of standard techniques in the failure analysis 
of the component, the treatment of the problem from the viewpoint of fracture 
mechanics. In using this approach the initial flaws or certain types of im- 
perfections which may exist in the material are treated as "cracks". These 
cracks may propagate under the fluctuating external loads as fatigue cracks, 
or under the sustained loads and sufficiently adverse environmental condi- 
tions as stress corrosion cracks, or a combination of both, namely as fatigue 
corrosion cracks. Aside from containing flaws and local damages in the form 
of notches, other stress raisers, or surface scratches which may be due to 
design, manufacture, or accidental causes, since pipes and containers are 
generally welded components they would almost always contain some imperfec- 
tions in the form of lack of total fusion, elongated slag inclusions, voids, 
arc burns, and a variety of other surface irregularities which may act as 
fracture nuclei. Thus, it appears that in actual structural components, the 
existence of defects which may lead to fracture initiation and propagation 
is unavoidable and therefore, one must design around them, primarily through 
the appropriate selection of material and dimensions of the component and by 
setting up rational defect inspection and acceptance standards. 

Fracture mechanics approach has been highly successful in dealing with 
this so-called fracture control problem which involves subcritical propaga- 
tion of fatigue and stress corrosion cracks and unstable fracture. The most 
widely used correlation parameter in this type of problems is the stress in- 
tensity factor which has been used quite effectively to predict the crack 
propagation rate in the component and for a given crack size, its residual 
strength. Using, then, this information, the "working life" of the component 
may easily be estimated. 

Up to now, fracture mechanics research as applied to pipeline and con- 
tainer safety, has dealt almost exclusively with the longitudinal flaws in 
the component. There is, of course, a good reason for this, namely that in 
pressurized cylinders under normal operating conditions, the hoop stress is 
by far the most dominant stress component. On the other hand, any secondary 
loads caused by known and unforseen factors would be primarily in the axial 
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direction. These are mostly bending and axial constraint stresses. Some 
of the sources of these stresses are support misalignment, variety of 
thermal fluctuations, ground settlement and earthquake, hydroelastic or 
aeroel Stic loading, axial constraint in buried pipes, gross bending in 
offshore piping caused by variety of oceanic disturbances and that gener- 
ated during installation, stresses arising from vehicular vibrations and 
impact, and the residual stresses inandaround the girth welds. One may 
also point out that the axial stresses would be much more time-varying than 
basically pressure-induced hoop stresses and hence would tend to facilitate 
the formation and propagation of fatigue cracks in the circumferential di- 
rection. If the thickness-to-radius ratio in the cylinder is sufficiently 
small, then a bending theory of shells may be used to formulate the cir- 
cumferential crack problem and, by approximating the internal stresses 
(such as residual and thermal stresses) by a linear function in the thick- 
ness coordinate, the general solution may be obtained by appropriately 
superimposing the results of membrane and bending solutions. The existing 
solutions of this problem take into account the membrane loading only and 
are based on an 8th order shallow shell theory in which the effect of the 
transverse shear stress is ignored and the crack surface boundary condi- 
tions are approximated by using the Kirchhoff assumption^ ^ [1-3]. The 
consequence of this assumption in shells is the same as in the bending 

' v 

of cracked flat plates [4-7]. It has been shown that if a plate bending 
theory compatible with the number of independent boundary conditions is 
not used to solve the crack problem, the computed stress Intensity factors 
may contain errors and the angular distribution of membrane and bending 
stress resultants would be different L4-7]. This would then make it dif- 
ficult to justify the superposition of membrane and bending solutions in 

7*1 . , . 

v That is, on the traction-free crack surface instead of satisfying 

the conditions Mp S = 0, V = 0 separately, the condition V+3M n r/9s = 0 is 
satisfied, thereby making the number of boundary conditions (four) and 
the order of the reduced differential equation (eighth) compatible. 

Here n and s are the coordinates normal and parallel to the crack, M ns 
is the twisting moment, and V is the transverse shear. 
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fatigue crack propagation and fracture analysis. 

The aim of this paper is to solve the problem of circumferentially 
cracked cylindrical shell by taking into account the effect of transverse 
shear and to obtain the stress intensity factors for bending moment as 
well as the membrane force as the external load. The problem for longitu- 
dinal crack was recently solved [8] by using a higher order shallow shell 
thoery described in [9]. In the present paper the effect of Poisson's 
ratio on the stress intensity factors and the nature of the out-of-plane 
displacement along the edges of the crack (i.e., bulging) will also be 
studied. 


2. FORMULATION OF THE PROBLEM 

In this problem the general shallow shell equations developed in [9] 
for an isotropic medium will be used. However, as shown in [10] and [11] 
for the classical shell theory, in this case too it can be shown that the 
related differential equations of the orthotropic shells car be factorized 
and reduced to the equations of isotropic shells if one assumes the material 
to be "specially orthotropic". 

Referring to Figure 1 and following [9] the equilibrium equations for 
a shallow shell may be expressed as 



(2.1) 





q(X r X 2 ) = 


( 2 . 2 ) 


H,i j • Vj 1 0 , (1 = 1,2, j = 1,2) 

' J * J 1 


(2.3) 


where N... M,,, and V., (1,j = l,2) are, respectively, membrane, moment and 

' J J 

transverse shear resultants, and the indicial notation and the summation 
convention are used. The components of strains are given by 


ij 


= i[U. . + U . . + Z .W . + Z .W .] , (l.j-1.2) 

£ ■ t J J > 1 » 1 f J 1 J » 1 


(2.4) 


where U-j , U£ and W are, respectively, X-j , X£» and Z-components of the 
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displacement vector, and the function Z^^) giving the equation of the 
middle surface is known. Let B-j and $2 be the angles of rotation of the 
normal to the shell surface. The (transverse) shear strains may then be 
expressed as 

9 i = w j + . 1 = 1.2 • (2.5) 
With the Hooke's law 


e 


ij 


a i jk^kt 


/h 


(2.6) 


and the relations giving M.j and in terms of 6^ and 0 ^, respectively, 
the formulation of the problem would be complete, technically the ten equa- 
tions (2.1) -(2.5) accounting for the ten variables e — , B^ » 9^ U^, and W. 

First, eliminating l)-| and U 2 from (2.4) one obtains the following 
compatibility equation: 


e ik e j£^ e i j,kt + Z ,ij W ,k«^ 0 


(2.7) 


(*) 


where e^ is the permutation symbol v . Next, defining the stress function 
F(X 1 ,X 2 ) by 


N ij e ik e j«. F ,k£ 


( 2 . 8 ) 


it is seen that (2.1) is satisfied and (2.2) and (2.7) reduce to 


+ Z ,ij e ikV.kt + q = 0 


e 1m e jn e kp e iq a ijkt F ,mnpq + hZ ,i j e ik e jfc W ,k£ 0 


(2.9) 

(2.10) 


Even for simple shell geometrip. for anisotropic materials the 


T*7 


e ll = 0 = e 22 * e 12 = 1 ’ e 21 * 


differential equations are not tractable. However, as in [10] if one assumes 
a special orthotropy, the related differential operators in these equations 
can be factorized and the problem can be made analytically tractable. Let 
the material be orthotropic and the stress strain relations be 


e ll ' hE-j ^ N ll " v 1 N 22^ 
e 22 = hl^ ( N 22 ‘ v 2 N 1 1 ^ 


_ 1 2 

e l 2 " IhG^ ’ 
E, " E„ 


( 2 . 11 ) 


In this case too the following is the condition for the factorization 
of the operators: 



1 + /v-j"v2 


( 2 . 12 ) 


The material satisfying the condition (2.12) is said to be "specially ortho- 
tropic". It has been shown that for certain orthotropic materials the valje 
of the shear modulus calculated from (2.12) and that measured experimentally 
are very nearly the same and consequently with the assumption of special 
orthotropy, the analysis can be simplified quite considerably [11]. If we 
now define 


E = /E-| , v = iA>i v£ , c = (Ei 


1 2 


(2.13) 


the stress-strain relations (2.11) become 

Ni 


J_ / II 

'11 hE 


1+v 


Ell _ ur (“72 ^ 22 ) » ^E » 


e 22 hE ( c2N 22"' ,N ll* 


(2.14) 


From the linear thickness variation of the stress components one obtains 

Mu = D (c 20 l,l + vB 2,2^ * M 12 = 2" ^ ^ B 1 ,2 + 8 2,1 ^ ’ 


Eh 3 


M 22 = D ^ ve l,l + 8 2,2 /c2) ’ D = 12(1 -v z ) 


(2.15) 
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Also, assuming a linear transverse shear stress-strain relationship, one 
finds 


9 1 " chB V 1 


0 = — v 
°2 hB V 2 


(2.16) 


where B is the effective shear modulus. Referring to [12], it will be 
assumed that 


5 _ E 

6 27TkT 


(2.17) 


Defining now the operator 

c c 3X^ C3X^ * 

equations (2.9), (2.10), and (2.3) may be reduced to 


v*v2F + hE( 
c c 

DV^W - (1 
c c 


3 2 Z 3 2 ? 3 2 Z 3 2 3 2 Z 3 2 v w _ n 

3Xf 3Xj ‘ 3X 1 3X 2 3X 1 3X 2 3X| TxJ' 

JLn2w3 2 Z_3l o ^ 2 Z 3 2 

Bh *c M axf " 3 x^X 2 3 x^X 2 


= 0 '^ 7 ‘ )q 


1+v 3 


36 o 96i 


R + -£!!L = r v 2g + (_£. 

P 1 3X 1 hB L Vl 2c 3X 2 ^3X 1 3X 2 ‘ 

„ t ii. i rv 2 8 +c i iv _l. ( !!i . !!i' 
e 2 3X 2 hB L c P 2 2 3X 1 V 3X 2 3X ] ' 


(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 
( 2 . 22 ) 


Equations (2.19) - (2.22) provide the formulation for an arbitrary shallow 
shell in terms of the unknown functions F, W, 6-j , and 02- 

Now let us assume that in the domain of interest the curvatures of the 
shell are constant. Then in (2.19) and (2.20) the terms involving Z may be 
replaced by 


3 2 Z _ J_ 3 2 Z _ _]_ 

3Xy " " R 1 * 3Xj ‘ ‘ R 2 


3 2 Z 


3X 1 3X 2 



(2.23) 
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Also, following [8], if one introduces the dimensionless quantities given 
in Appendix A, equations (2.19) - (2.22) may further be simplified as 


v “* • X 7 (X 1 ^ ' 

2X 12 3x3y + X 2 ■ 0 

(2.24) 

V 4 W + X 2 ( 1 -<V 2 ) ( 

& - 2X.2 32 t , * 

3 y 2 12 3x3y 2 3x z/ 



•X-(l-«72)iq , 

(2.25) 

(1-^) 6x+ |^= K 

1+v 3 

2 3y 3x " 3y ’ 

(2.26) 

(1-kV 2 )8 + = k 

v ' y 3y 

,. . 38 38 

1+v 3 , p x _ 

2 3x v 3y 3x 

(2.27) 


The constant a used in Appendix A to normalize various quantities is a charac- 
teristic length parameter in the shell. Usually in crack problems the shell 
is assumed to be "infinitely large" and a is taken to be the half crack 
length. 

Defining now a new function n by 

38 38 

- TT - 1* (2 ' 28) 


equations (2.26) and (2.27) may be expressed as 

a = r o2a _ _ 1+v 3 n 

0 x 6 x 3x k 2 3y * 


a = <7 2o . iSt + t hX M 

^ 6 y 3y K 2 3X • 


and, it may easily be shown that 


V 2 8 = — + — -P + J - 

x 3y 3 x z 3x3y 


3 2 8„ a 2 6. 


V z 8.. = 


30 


3 Z 8„ 3^8 


3a 3x3y 3y 




(2.29) 


(2.30) 
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Also, if we define the function ip by 

*(x,y) • + - 5 ^) - w , (2.31) 

from (2.29) and (2.30) we obtain 

(2.32) 

In (2.32) eliminating a and then using (2.31) we find 

<V 2 ^ - ip - w = 0 . (2.33) 

Similarly, eliminating w, (2.29) yields 

k ^ v 2 n - n = 0 . (2.34) 

The solution of the shell problem must then satisfy the differential equa- 
tions (2.24), (2.25), (2.33) and (2.34) and all the necessary boundary con- 
ditions. 

3. CYLINDRICAL SHELL WITH A CIRCUMFERENTIAL CRACK 

Consider now a cylindrical shell containing a circumferential through 
crack of length 2a shown in Figure 2. In this case X-j = 0 = X^, and if 
we further assume that q = 0, (2.24) and (2.25) become 

v 4 <j> - (X,/X ) 2 “7 = 0 , 

v-w + (XX 2 )2(1 -kV2) 0*0. 

Eliminating <t>, from (3.1) and (3.2) it follows that 

v 4 v 4 w + x£(l-<v 2 ) I^Sr = 0 . (3.3) 


(3.1) 

(3.2) 


6 = M + K hy. 22 

ax 2 ay 

6 = 3i l-v M 

y ay 2 ax 
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The problem will be solved by using Fourier transforms. It will be assumed 
that through a proper superposition the original shell problem has been re- 
duced to a perturbation problem in which self-equilibrating force and moment 
' resultants acting on the crack surfaces are the only nonzero external loads 
(hence, the assumption q = 0). Thus, in some neighborhood of the crack in 
“which the stresses are expected to be nonzero, the transform of w and its 
inversion may be expressed as 


f(x,a) 


w(x,y)e lniy dy 


(3.4) 


w(x,y) 


2tt 


f(x,a)e uy da 


(3.5) 


Substituting from (3.5) into (3.3) one obtains an 8th order ordinary linear 
differential equation in f. Looking for a solution in the form f=R(a)exp 
(mx), the characteristic equation of the problem is found to be 

m 8 - (k X 2 + 4a 2 )m 6 + (6a 4 + x£ <X^a 2 )m 4 


- 4a & m 2 + a 8 = 0 


(3.6) 


If we define 


p = m 2 - a 2 (3.7) 

equation (3.6) can be written in terms of p as follows: 

p 4 - icX^p 3 + X^(l-2<a 2 )p 2 + Xij(2a 2 - <a 4 )p + X^a 4 = 0 (3.8) 

It should be emphasized that the roots of (3.8) are in general complex and, 
of course, are not known as a function of a in closed form. Appendix B 
describes a convenient procedure for solving (3.8). After solving (3.8) 
let the roots of (3.6) be ordered such that 

Ro(mj<0 , nij +4 = -m^ , j = 1 , . . ,4 (3.9) 
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The solution f(x,a) of the resulting, ordinary differential equation sat- 
isfying the regularity conditions at x = + « may then be expressed as 


r 4 


f(x,a) = < 


l R • (a)exp(m .x) , x>0 , 

j=l J J 


t R.(a)exp(m.x) , x<0 . 

U-5 J J 


If we let 


*(x,y) = ^ 


g(x,a)e" 1ya da , 


from (3.1), (3.2) and (3.10) we find 


g(x,a) = < 


(X 2 /X)2 V (mj/Pj) 2 Rj( a )exp(mjX) , x>0 , 


(x 2 /x) 2 l (mj/Pj) 2 Rj( a )exp(mjX) , x<0 


Similarly, assuming 


n(x,y) = 27 


h(x,a)e iya da 


from (2.34) we obtain 

* 

A-j (a)exp(r 1 x) , x>0 

h(x,a) = < 

A 2 (a)exp(r 2 x) , x<0 

where 


r l " • r 2 = ^° 2 + kO-.]' 1 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


Also, let us assume that 


e(x,a)e 1ya da . 


(3.16) 


/ 


'l'(x.y) 


2 TT 


—co 


It can be shown that the remaining differential equation ( 2 ; 33) is satisfied 
if e is assumed to be 


4 R (o) 

l „ — t exp(m.x) 
i kPj-i J 

8 R.(a) 

l sbr exp(m j x) 

l. 3 J 


x>0 

x<0 


(3.17) 


The expressions given by (3.10), (3.12), (3.14), and (3.17) satisfy 

t*\ 

the differential equations of the problem' '. If one now determines the 
arbitrary functions Rj(a), (j = l,..,8) and A-j(a), (i = 1 ,2) in such a way 
that the boundary conditions of the problem are i.:so satisfied, one then 
has the solution. 


4. BOUNDARY CONDITIONS 

As mentioned earlier, the only external loads in the problem are the 
self-equilibrating force and moment resultants on the crack surfaces. These 
forces can be decomposed in such a way that in solving the problem one needs 
to consider (with respect to x = 0 plane) either symmetric or antisymmetric 
loading only. In these two cases the following conditions of symmetry will 
be satisfied: 


N xx (x,y) = N xx (-x,y) , N xy (x,y) = -N xy (-x,y) , 

M xx ( x *y) = M X x ( - x ’ y) , M xy (x,y) = -M xy (-x,y) , (4.1) 

V x (x,y) = -V x (-x,y) , 

v Mt should perhaps be pointed out that (3.17) is a solution satisfying 
the differential equation (2.33), and is not the most general solution. 
However, with (3.17), since the solution thus found satisfies all the dif- 
ferential equations and, as will be shown, all the boundary conditions, 
it must be the solution of the physical problem. 
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for the symmetric problem, and 

N xx ( x »y) = -N xx (-x,y) , N xy (x,y) = N xy(" x,y ) * 

M xx ( x > y ) = -M xx (- X » y ) • M X y( x ’ y ) = M xy (-x»y) » (4.2) 

v x U,y) = v x (~ x »y) > 

for the antisymmetric problem. One may note that (4.1) and (4.2) are valid 
for all values of x and y, and in the odd functions the discontinuity at 
x = 0 may be allowed only on the crack surface, outside the crack all these 
functions (indeed, all field quantities) must be continuous. It is there- 
fore clear that in solving the problem one needs to consider only one half 
of the cylinder, say x > 0. Symmetry conditions similar to (4.1) and (4.2) 
are, of course valid for all the remaining field quantities. Thus, there 
are only five unknown functions R^,..,R^, and A-j which may be determined 
from five conditions specified at x = + 0. 

Consider now the symmetric problem *or a circumferentially cracked 
shell. Noting that before the superposition which led to the perturbation 
problem the crack surfaces were free from all external loads, and since 
outside the crack all quantities are continuous, from (4.1) it may be con- 
cluded that 


N xy (0,y) = 0 , M xy (0,y) = 0 , V x (0,y) = 0 , -»<y<» , (4.3) 


and N (+0,y) and M (+0,y) are known functions in -1 <y< 1 . The problem is 

AA AA 

then a mixed boundary value problem. Using (4.3) three of the five unknowns 
(Rl , . . ,R^,A-j ) can be eliminated. The remaining two may be determined either 
from a system of dual integral equations or a system of singular integral 
equations arising from the mixed boundary conditions. 

In terms of the normalized quantities the relevant force and moment 
resultants are given by 



3 2 4 > 

3x9y ’ 
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M. 


xx 

M xy 

3V., 


_a 

hX 1 


33 33 

(_* + v A) 

' ^ v/ Jy / 


3X 




1 96 

a IziL ( JjL + - 11 ) 

hx^ 2 ' 3y 3x ' 

3 2 33 

W ^ X 


3y 3x3y 3y 


(4.4) 


Using (4.4) and the results obtained in Sections 2 and 3, for x>0 these 
quantities may be expressed in terms of R-|,..,R 4 and A-j as follows: 

m-x 


1 X« 2 ** III' 2 111.5A ■; 

N xx (x,y) = - 2? [(->-> “ 2 I R j<“> e ' d “ 

_ 0 O 1 J 

. r° X 9 2 4 m-X • , 

N xy (x,y) = 27 C(T> ^ ^ Rj(a)e J ]e _lya da , 

— oo ^ J 


"»<*•»> ■ h hT 


K ( 1 ~ v ) 


2 , 2 

c 

00 4 m.-va 2 


V„- -iyc 


I ~ j, T R i( a ) e J e 
f KPj-1 J 


da 


2 r° 


r l X -iya. 


iar^A^ (a)e e y da 
4 m.R.(a) m.x 


l ad r° 4 m. 

io H^r eJ 

— OO 1 J 

["(a 2 + rf )A, (a)e 


e- ly “da 


1 a< 
2tt hx 4 


r ' V 1y ° 


da 


— V 

3y 


K P ,m. 


^V y “da 


X ( x,y ) = - 27 j ia \ ^f R j (a)e 

“oo ^ 

1 1 -v f 2 , n r l X -iya . 

- 27 tc ~ 2 ~ a 2 A^ (a)e e J da 


(4.5) 

(4.6) 


(4.7) 


(4.8) 


(4.9) 


We will attempt to solve the problem by reducing it to a system of 
singular integral equations. The problem is "symmetric" and the "normal" 
membrane and moment resultants are specified on the crack surface. There- 
fore, the natural "dual" quantities which should be considered as the new 
unknown functions are 
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(4.10) 


~ u(+0.y) - G,(.) , i 6 x (+0,y) = G 2 (y) , 

corresponding to the "normal" displacement and rotation on the crack surface, 
respectively. In (4.10) the derivatives of the dual quantities are used to 
make them dimensionally consistent with N and M vv and in order to insure 
that the resulting integral equations will have Cauchy type singularities 
[13]. The mixed boundary condition along x = 0, -~<y«» may now be expressed 
as follows: 


1 1m N (x,y) = F, (y ) , 

x++0 xx 1 

-l<y<l . 

(4.11) 

lim M (x,y) = F ? (y) , 

x-h- 0 xx 

-1 <y<l , 

(4.12) 

u(0,y) = 0 , l<|y|<» , 


(4.13) 

6 x (0,y) = 0 , 1 < |y l <°° 

• 

(4.14) 


Referring to the definitions (4.10) and the physical conditions (4.13) and 
(4.14), it is clear that the functions G-| and G2 must satisfy the following 
single-val uedness conditions: 


rl 

G^yjdy = 0 

-1 


G 2 (y)dy = 0 . 

-1 


From (2.4), (2.14), (2.8), and (2.1) it may be shown that 


(4.15) 


-~z u(+0.y) ■ - *(+0,y) + (x 2 /a) 2 ^ w(+0,y) . (4.16) 

Then, after some manipulations, the quantities defined by (4.10) may be ob- 
tained as follows: 


G 1 (y) 

G 2 (y) 


lL ) 2 j_ r 1 

' A 2w 


, 14 m? 

a \ <V pf )R J (o)e ' 

J 


■lya 


da 


1 

2tt 


f" 4 rn.R.(a) 


« l 


_J_ A. 


■V 


e-' ya do 


(4.17) 


00 
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(4.18) 


2tt 


T~ ( 

*-oo 


2 A-j (a)e 1 ^ ot da 


Now, using (4.3) and inverting (4.6), (4.8), (4.9), (4.17), and (4.18), 
we obtain five linear algebraic equations in R^,..,R^, and A 1 which, after 
some modifications and by using (4.13) and (4.14), may be expressed as 



(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 


The solution of the system of equations (4.20) - (4.23) is given in Appendix 
C. By substituting from (4.19), Appendix C, (4.5) and (4.7) into (4.11) and 
(4.12), we obtain two integral equations to determine the unknown functions 
G-| and G 2 of the following form: 


,1 

1 im 

X-++0 i - ) 


l G;(t)dt [ H..(a,x)e 1(t ' y)a dc. « F. (y) 
1 J J 1 J ' 


» 


i = 1,2, -1 <y<l (4.24) 

By examining the functions H.., it can be seen that they are bounded every- 

1 J 

where in -~<a«». Therefore any possible singularity of the kernels in (4.24) 
at y=t must be due to the behavior of H-jj(a.x) as a-++°°. Note also that H^. 
contains exponential damping terms of the form exp(mjx) and exp(r^x), where 
Re(mj)<0, Re(r-j)<0. However, since in limit x will go to zero, for y= t 
this damping does not insure the convergence of the inner integrals in (4.24). 
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The major difficulty in this problem, of course, is that the functions 
m.(a) are not known explicitly in terms of a. For the purpose of examining 

J 

the singular behavior of the kernels in (4.24) and for extracting the 
singular parts, all one needs, however, is the asymptotic behavior of m. 

s) 

and r-j as |al-*». Thus, from (3.6) - (3.8) and (3.15) it can be shown 
that for large values of |a| we have 


P • P ■ 

m j(o) ■ -M (1 + i£z - gjr + 

r l ( “> = -l«l o + K (il v ) 0 ^ - • 


) , 

(4.25) 

) . 

(4.26) 


Using now the relations (4.25) and (4.26) and separating the asymptotic 
values of for large |a|, the kernels in (4.24) may be expressed as 


|V (t ' ,)8 - 


* f H”.(a,x)e^ t " y ^ a da 


10 


- H^ j (a.x)]e 1(t ‘ y, “<fa 


(4.27) 


where H*.. is the asymptotic value of for |a|-*». On the right hand 
side of (4.27) the first term gives Cauchy type kernels 1 / ( t-y ) on the main 
diagonal terms, and the second integrals are uniformly convergent for all 
t and y (in which, the limit x = 0 can therefore be put under the integral 
sign). After the asymptotic analysis and some lengthy but straightforward 
manipulations the integral equations and the kernels may then be expressed 
as follows: 


rlG,(t) * i 

j ■ dt + l k-j j(y,t)Gj(t)dt = 2 1 iF 1 (y) , 

2 Go(t) 2 rl 
~ 7 T [ T7* + l J k 2j<y.t)Gj(t)dt 


-1 <y<l , (4.28) 


= 2ir a F 2 (y) . -i<y<i 


(4.29) 
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(4.30) 

(4.31) 

(4.32) 

(4.33) 

where the functions N.(a) and Q.(a) are giver in the Appendix C. Using the 

J J 

results of Appendix C it can be shown that k^(y»t) = ^(y.t). 

5. THE ASYMPTOTIC STRESS FIELD AROUND THE CRACK TIPS 

The solutions of the singular integral equations (4.28) and (4.29) have 
integrable singularities at the end points y = +1 and are of the following 
form: 


r» 2 4 2m? 

k 11 (y,t) = J [(y) a 2 y. -pf Qj(a)-l]sina(t-y)da , 


k 12 (y,t) 
k 21 (y,t) 
k 22 (y,t) 


x„ 2 f o° 4 m? 


A 9 r mi 

2(y) a 2 l -4 N .(a)sina(t-y)da , 

A J 0 1 p j J 

2x 9 r 4 i 

- TIT y i (m? - va 2 )m‘tO i (a)sina(t-y)da , 
A J 0 1 Pj J J J 

00 4 X** 

-2 y j [I FJ (m? - va 2 )m'jN j (ct) 


1-v 2 


<(1 -v) 2 ar^ + — 2 — ]sina(t-y)da 


G^y) = g^yJO-y 2 )’ 1 * » i = i»2 , (5.1) 

where g-j and g 2 are bounded in -l<y<l. Similar to the plane problems, it 
can be shown that the behavior of the stress distribution in the immediate 
neighborhood of the crack tips is dependent on g.(+l) only. To show this, 
one needs to substitute (5.1) through the expressions of R-|,...,R 4 given 
in Appendix C and A-| given by (4.19) into the original expressions for 
and M.., (i,j = x,y), such as those given by (4.5) - (4.8). Using now the 

* J 

relation [14] 


u-1 -ba r sin 
a e 


{ cos ,(ca)da = - (b2 ™,* <cos )(,jtan " 1 5 ) 


(b>0 ,y>0) , (5.2) 
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and [8, 10,11] for large values of |a 


f jiiL 

-1 /TF 


+ 


e lat dt = { g( 1 )exp[i (a - J sign(a))] 

g( -1 )exp[-i (a - J- sign(a) )] + 0 (-^y)} , 


(5.3) 


around the end point y=l, x = 0 the leading terms of the asymptotic stress 
and moment resultants may be expressed as 


N xx (x ’ y > 

V*,y) 

N X y( X »y) 

M xx (x > y) 

M yy (x,y) 

M xy U.y) 


f° — (1 + a| x| )e~ ' X ^ a sin[(l-y)a - J]da 
f — (1 - a|x| )e’^ a sin[(l-y)a - J]da 

* n 


2/27 

— 0 + 
0 /a 

9,0) 

2^7 J 

O 8 

9,0) 

2/27 

.OO 

/a xe" 

0 

g 2 0) 

h 

r _L 

2/27 12a J 

0 G 

g 2 0) 

h 

_L 

2/27 12a 

0 /a 

g 2 0 ) 

_h_ 

o _ 



— (l + a|x|)e'l x|o sin[(l-y)c. - j ]da 
/a 

f J- (1 -a|x|)e' |x Kln[(l-y)« - J]da 

J i/nT 


2^27 


(5.4) 

(5.5) 

(5.6) 

9 

(5.7) 

9 

(5.8) 

(5.9) 


From (5.4) - (5.9) it is seen that, aside from {he magnitudes as rep- 
resented by g-j(l) and g^O), the asymptotic behavior of the membrane and 
bending stress distributions around the crack tip will be identical. This 
is, of course, in agreement with the uncoupled in-plane and bending results 
for flat plates [4-7]. Defining the polar coordinates r,e by 


x = rsine , y - 1 = rcose , 


(5.10) 


evaluating the integrals, and observing that in dimensionless quantities 
(see Appendix A) 
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. (I.j * x,y) , 


(5.11) 




12 az M 

l M* • — o • • 
h i j i j 


from (5.4) - (5.9) the asymptotic stress distribution may be obtained as 


° xx (r,e >Z ) = 


gi(D + zg 2 (l) . 8 , 


2 / 2 r 


r j o i uon 

[4 cos 2 - 4 cos y] , 


(5.12) 


9i ( ^ ( 1 ) o a 1 c Q 

O (r,e,z) = — [| cos | + i cos y] , (5.13) 

yy 2^2r * L * 6 


o xv (r,e,z) . - 9 - 1 - ^ 01 [- } sin f ♦ ' sin . (5.14) 

xy 2 /2r d 4 


In this "symmetric" problem the stress intensity factor is defined by 


k-| = lim /2( X 2 -a) a-|i(0,X 2 ) 


(5.15) 


which, referring to Appendix A and (5.12) may be obtained in terms of g-j(l) 
and g 2 (l ) as follows: 


k,(Z) - - [g,(l) + f 9 2 0)] • 


(5.16) 


It should be noted that even though the analysis given in the previous 
sections is valid for specially orthotropic as well as isotropic shells, 

(a) because of the dependence of x and y on the orthotropy constant 
c = (E-|/E 2 )\ r and 0 defined by (5.10) are not the physical coordinates, 
and consequently the angular distribution of the stress state in orthotropic 
shells would be different from that of isotropic shells as given by (5.12) - 
(5.14); and (b) since the roots and r^ are very heavily dependent on c, 
the numerical results obtained for the isotropic shells by assuming c = 1 
cannot be eadily adopted to the orthotropic shells. 


6 . CRAC.. SURFACE DISPLACEMENTS 


After obtaining the functions G-| and G 2 upon solving the integral 
equations (4.28) and (4.29), through the expressions for R-|,..,R^, and 


I 


A.| , any desired field quantity in tre shell can be expressed in terms of 
G-j and G 2 and be easily evaluated. One such group of quantities of some 
practical interest is the components of crack surface displacement vector. 
In the symnetric problem under consideration, referring to the definition 
(4.10), the in-plane component of the crack surface displacement (i.e., 
the crack opening displacement) may easily be obtained as 


u(+0,y) 


f l 

G.|(t)dt , - 1 <y< 1 

y 


(6.1) 


The component of the displacement vector which is perpendicular to the 
shell surface, i.e., w(x,y), is given by (3.5) and (3.10). Expressing 
again R-|,..,R^ in terms of G-j and G 2 (Appendix C), w may be obtained as 
follows: 


w(+0,y) 



G-j (t)dt 


.00 4 

l 

0 1 


Q 4 (a)sinot( t-y)da 


G 2 (t)dt 


.00 4 

l N .(a)sina(t-y)da 
0 1 J 


» 


-1 <y<l 


(6.2) 


where Q. and N- are given in Appendix C. 

J 


7. NUMERICAL SOLUTION AND RESULTS 

The singular integral equations (4.28) and (4.29) subject to single- 
valuedness conditions (4.15) are solved by using Gauss-Chebyshev integra- 
tion formulas. Thus, equations (4.28), (4.29) and (4.15) are, respectively, 
replaced by 



LWvVSm'V 1 * 2, W 

m= I 


i = 1 , . . ,n-l 


(7.1) 
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(7.2) 


n , ; 

l V'TT 

0=1 0 x 


‘o- y > 


m= 1 


= 2v J F 2 (y i ) , i = 1 . ,n-l , 


kW )s ° ’ J I 1 W J g 2 (t J ) 


= 0 , 


where 


tj = cos (jrf ") » j = 1 »n , 


-2i-l 


y i ~ cos( 2 n _2 ^ f i - 1 , . . ,n-l * 


W 1 * M n = 2Tn'-"lT ' W j = n-1 ’ j' 2 ---*"- 1 


(7.3) 


(7.4) 

(7.5) 

(7.6) 


As a numerical example a cracked cylindrical shell under uniform mem- 
brane and bonding loads is considered. To make the practical applications 
of the results more convenient, these two loads are considered separately. 
The calculated results for the stress intensity factors are normalized with 
respect to the corresponding flat plate values. For example, if the applied 
loads are 


N n (0,X 2 ) = -N n = -ho m , M 11 (0,X 2 ) = 0 , -a<X 2 <a , (7.7) 

the input functions in the integral equations (4.28) and (4.29) become 


F-j (y) - -N n /hE , F 2 (y) = 0 


(7.8) 


and the corresponding flat plate stress intensity factor is a m /a\ Then the 
membrane and bending stress intensity ratios k^ and k^ m are obtained from 
(5.16) as follows: 



(7.9) 
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k 1 (h/2)-k 1 (0) E 


o m i/a 
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2 o 


m 


21 g 2 (1) 


(7.10) 


The "bending stress intensity factor" is thus based on the bending stresses 
at the outer layer of the shell, Z = +h/2. 

Similarly, if the external loads are 

N n (0,X 2 )=0 , M 11 (0,X 2 ) = -M n = - ^ o b , -a<X 2 <a (7.11) 

Then the Input functions become 

M n 

F-|(y) = 0 , F 2 (y) = - ^ , -1 <y<l . (7.12) 


In this case too the corresponding flat plate stress intensity factor is 
defined by a^/a, and the stress intensity factor ratios are obtained as 
follows: 


M°) r 

k mb ‘ = ' V 1 ' • 

_ k l th/2) - k,(0) _ E h m 
c b ^a ‘ ‘ 2a b 2a 9 2 (1) 


(7.13) 

(7.14) 


In the numerical calculations the effective transverse shear modulus for 
the shell is assumed to be B = 5G/6, G being the shear modulus of the material. 
Also, the Poisson's ratio is taken as v=l/3 in all calculations except for 
one set of results where the effect of v is investigated. One may note that 
*2 is used in this analysis is the standard shell parameter A defined in the 
formulation of the problem by using the classical (i.e., the 8th order) shell 
theory. Also note that in the present analysis there is an additional para- 
meter a/h, which, within the confines of the shallow shell theory, gives the 
thickness effect. 

The numerical results are shown in Tables 1-5. To help visualizing 
the trends some of the results are also shown in Figures 3-6. Tables 1-4 
show the stress intensity factor ratios k^, k bm , k mt) , and k bb defined by 
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(7.9), (7.10), (7.13) and (7.14), respect vl ly . In these tables the 
values given for A 0 = [12(l-v 2 )] l< a//Rh^ = 0 correspond to the flat 
plate. Both *£ and a/h have been varied from 0 to 10. For >10 and 
a/h > 10 the linearized shallow shell theory used in this study is prob- 
ably not valid. For very thin shells (i.e., a/h = 10) under uniform 
membrane loading (Tables 1, 2 and Figures 3, 4) it is found that the 
membrane component of the stress intensity factor k nm is indistinguish- 
able from the results found by using the classical theory [1]. However 
the bending stress intensity factor k bm is quite different. The results 
given in Tables 1-4 show that, particularly for large values of Ao, the 
thickness parameter a/h may have considerable effect on the stress in- 
tensity factors. Table 3 and Figure 5 show that when A^-^O the results 
are in good agreement with the flat plate bending results given in [5] 
and [6] and the axially cracked shell results given in [8]. 

For a specific geometry a^ = 4 and a/h = 2, Table 5 shows the effect 
of the Poisson's ratio v on the stress intensity factors. For the two 
most important components k^ and k bb (i.e., for the primary stress in- 
tensity factors under membrane and bending loads) the effect of v does 
not seem to be significant. Hence the results given in Tables 1-4 and 
calculated for v=l/3 can be used for materials with a Poisson's ratio 
0.2 <v <0.4, which may cover nearly all structural materials. 

Figures 7 and 8 show some sample results for the out of plane dis- 
placement w(+0,y) calculated on the neural surface along the line of 
the crack. The trend seems to be quite similar to the bulging results 
obtained for the axially cracked shell [11]. It should again be em- 
phasized that in these numerical calculations the crack surface membrane 
stress is compressive and the bending moment is applied in such a way 
that the outer layer of the shell on the crack surface is again compres- 
sive. Under these loads, as seen, from the figures, on the crack surface 
the displacement w is in the outward direction. 
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APPENDIX A 

The Dimensionless Quantities 
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APPENDIX B 

Procedure for Solving a Quartic Equation 

Consider the following quartic equation: 

p 4 + a 3 p 3 + a 2 p 2 + a-|P + a Q = 0 (B.l) 

Although in the shell problem the coefficients a Q ,..,a 3 are real, for the 
sake of generality, here it will be assumed that they are complex. There- 
fore the roots of (B.l) are in general complex. Assume that (B.l) can be 
written as the difference of two squares: 

(p 2 + Ap + B) 2 - (Cp + D) 2 = 0 , (B.2) 

Then, from (B.2) it follows that 

p 2 + (A + C)p + (B + D) = 0 , (B.3) 

p 2 + (A - C)p + (B - D) = 0 . (B.4) 

If the constants A, B, C, D can be determined in terms of aj, a 2 » a-| , a Q , 

then the four roots can be found in a straight-forward manner by solving 

equations (B.3) and (B.4). 

Comparing (B.l) and (B.2), one can write: 


2A = a 3 , 

(B.5) 

2B + A 2 - C 2 = a 2 , 

(B.6) 

2AB - 2CD = a ] , 

(B.7) 

B 2 - D 2 = a . 

0 

(B.8) 


Eliminating C and D in (B.6) - (B.8) and defining 2B - u, one finds: 
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(B. 9) 


u 3 + b,,u 2 + b,u + b„ = 0 
2 1 o 


where 


b 2 = -a 2 


b l ‘ 


a, a., - 4a„ 
13 o 


= 4a o a 2 


1 2 _ 


a. a: 


‘1 o 

o 3 


(B. 10) 


Let u-j be a root of the cubic equation (B.9), (see, for example [15] for 
determining the roots of a cubic). Then 



(B.ll) 


and from (B.5) 



(B.12) 


Once A and B are known, one can determine C and D from (B.6) and (B.8) as 
follows: 





a 2 


(B.l 3) 
(B. 14) 


It is seen that C and D are multiple-valued functions and one should choose 
the correct branch in order to obtain the correct solution. 

C 2 and D 2 can be written in the complex form as follows: 

ie, ie~ 

C 2 = P] e , D 2 = p 2 e c (B. 1 5) 


where p-j , p 2 , 0-| and e 2 can be determined by using (B.13) and (B.14). Thus 


ie,/2 

^Te ' , 

i ( e , / 2 + tt ) 

C 2 = /p] e 


i 0o/2 

i ( Gp/ 2 + tt) 


2 , 

D 2 = ^02 e 

(B.16) 
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All pairs of (C,D) satisfy equations (B.6) and (B.8) but not (B.7). 
Therefore C and D must be selected such that equation (B.7) is satis- 


fied, i.e.. 
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APPENDIX C 


Expressions for R.(a), j = 1 , . . ,4 
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The expressions of N^ and Q 2 are obtained from (C.2) and (C. 
the indices in m. and p. sequentially from 2, 3, 4 to 3, 4, 

J J 

for and Q 3 the indices are replaced by 4, 1, 2, and for N 
1, 2. 3. 


(C.l) 


(C.2) 


(C. 3) 


(C.4) 
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4 and Q 4 by 
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a/h = 0.5 

a/h = 1 

a/h = 2 

a/h = 5 
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Table 1 Stress intensity factor ratio k 

rrm 


Table 2 Stress intensity factor ratio k^ 


X 2 

a/h = 0.5 

a/h = 1 

a/h = 2 

LT> 

II 

.c 

a 3 

a/h = 10 

0.0 

- 

0.000 

0.000 

0.000 

0.000 

0.5 

- 

0.050 

0.047 

0.046 

0.046 

1.0 

0.101 

0.100 

0.098 

0.098 

0.099 

2.0 

0.103 

0.110 

0.119 

0.130 

0.139 

3.0 

0.066 

0.052 

0.052 

0.073 

0.089 

4.0 

0.036 

-0.010 

-0.037 

-0.024 

-0.004 

5.0 

0.017 

-0.055 

-0.117 

-0.126 

-0.106 

6.0 

0.002 

-0.087 

-0.179 

-0.217 

-0.202 

8.0 

-0.018 

i 

o 

u 

o 



-0.264 

-0.359 

-0.359 

10.0 

- 

-0.161 

-0.322 

-0.463 

-0.481 
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Table 3 Stress intensity factor ratio k. L 

bb 


x 2 

a/h = 0.5 

a/h = 1 

a/h = 2 

a/h = 5 

a/h = 10 

0.5 

- 

0.737 

0.693 

0.660 

0.653 

1.0 

0.738 

0.672 

0.636 

— 

0.612 

0.608 

2.0 

0.648 

0.549 

0.508 



0.493 

1 

0.493 

3.0 

0.615 

0.484 

0.422 

0.401 

0.401 

4.0 

0.602 

0.455 

0.376 

0.343 

0.342 

5.0 

0.595 

0.441 

0.350 

0.308 

0.303 

6.0 

0.592 

0.434 

0.335 

0.284 

0.277 

8.0 

0.589 

0.426 

0.319 


mm 

10.0 


0.423 

0.311 

0.236 

0.219 


Table 4 Stress intensity factor ratio k mb 


X 2 

a/h = 0.5 

- 

a/h = 1 

a/h = 2 

a/h = 5 

a/h = 10 

0.0 

- 

0.000 

0.000 

0.000 

0.000 

0.5 


0.015 

0.013 



1.0 

0.034 

0.032 

0.028 

0.025 


2.0 

0.044 

0.047 

0.045 

0.041 

0.039 

3.0 

0.341 

0.046 

0.046 

0.044 

0.042 

4.0 

0 036 

0.042 

0.043 

0.041 

0.040 

5.0 

0.032 

0.038 

0.039 

0.038 

0.036 

6.0 

0.029 

0.035 

0.036 

0.035 

0.033 

8.0 

0.024 

0.030 

0.031 

. . 

0.030 

0.029 

10.0 

- 

0.026 

0.028 

0.027 

0.026 
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Figure 7 Displacement component perpendicular to the shell surface, 
w(0,y) . N xx t 0, M xx = 0, a/h = 5, = 4 




